Spectral properties of high-gain 
parametric down-conversion 



K. Yu. Spasibko 1 , T. Sh. Iskhakov 2 , M. V. Chekhova 2 1 

1 Department of Physics, M. V. Lomonosov Moscow State University, Leninskie Gory, 119992 

Moscow, Russia 

2 Max Planck Institute for the Science of Light, Gunther-Scharowsky-Strafie 1/Bau 24, 91058 

Erlangen, Germany 
drquantum @ hotmail com 



Abstract: High-gain parametric down-conversion (PDC) is a source of 
bright squeezed vacuum, which is a macroscopic nonclassical state of light 
and a promising candidate for quantum information applications. Here we 
study its properties, such as the intensity spectral width and the spectral 
width of pairwise correlations. 
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Bright squeezed vacuum (BSV), the state generated via high-gain parametric down- 
conversion (PDC) and having a large number of photons per mode, is currently in the focus 
of attention in quantum optics and quantum information. The reason is that BSV is considered 
as a state that can manifest macroscopic entanglement Q. At the same time, because BSV can 
have huge numbers of photons per mode, it can be easily involved into nonlinear optical inter- 
actions |2] and interactions with material quantum objects such as atoms, molecules, quantum 
dots etc. In particular, its application to quantum memory is discussed [3]. Recently, BSV has 
been used for quantum imaging [4] and absolute calibration of photodetectors Q. Because 
of these applications, it is very important to know the spectral properties of BSV emitted via 
high-gain PDC. 

It should be noted here that the spectra of PDC are essentially two-dimensional (see, for 
instance, Ref. O): the shape of the frequency spectrum depends on the selected angle of emis- 
sion, and vice versa. In this connection, it is worth mentioning an interesting concept of two- 
dimensional coherence of PDC radiation (7) . Therefore, a general consideration should involve 
both frequency and angular spectra. However, in this work we only consider frequency spectra 
and all measurements are performed with an aperture selecting only collinear emission. 

At high gain, PDC does not any more generate two-photon light, as is the case of low gain. 
The radiation of high-gain PDC contains not only photon pairs, but four-photon states, six- 
photon states, and so on, and the mean photon number per mode is large. The mean photon 
number per mode is in one-to-one correspondence with brightness', this is why the state can 
be called bright squeezed vacuum. Under these conditions, correlations in SV cannot be ef- 
ficiently studied by measuring the number of coincidences. Instead, one should measure the 
noise (variance) of the intensity difference for the two selected beams l8l[9lfT0l. 

The spectral and angular properties of high-gain PDC (known earlier as parametric super- 
fluorescence full ) have been well studied theoretically EU [121 H3 . In particular, it has been 
known for a long time ifTTIl that the spectral width of both signal and idler beams should increase 
with the gain. Physically, this is clear from the fact that at high gain, there is exponential am- 
plification of down-converted light along the crystal. Therefore, effectively it is not the whole 
crystal that contributes into the signal and idler intensity but only its exit part, where the sig- 
nal/idler intensities (photon numbers) are most high. For the same reasons, the spectral width 
of correlations, usually determined by the inverse pump pulse width/coherence time, should 



also increase at high gain. Similar behavior should be observed for the angular spectra: both 
the total angular width and the far-field zone speckle size (correlation radius) should increase at 
high gain. Indirect experimental evidence for this was observed in Ref. [8J. The speckle size has 
been directly measured as a function of the parametric gain in Ref. iTPfll . However, no measure- 
ments on the frequency spectra of high-gain PDC have been performed, although there was an 
experiment on the study of spectral correlations for a similar effect, four- wave mixing fT5l . 

In this work, we study the spectral properties of BSV generated via high-gain PDC from a 
picosecond-pulse pump. Due to a rather short coherence time of the pump, we are able to re- 
solve the spectral width of correlations between signal and idler radiation. By scanning both the 
signal and idler frequencies separately, we can see the two-dimensional distribution of correla- 
tions, similar to the 'joint two-photon intensity' observed for two-photon light 021 021 OS • We 
also measure the spectral width versus the parametric gain and observe the spectral broadening 
in accordance with ifTTIl . Finally, we observe some interesting effects caused by super-bunching 
in high-gain PDC. 

The experimental setup is shown in Fig.Q] BSV was obtained in a 2 mm long BBO crystal 
cut for type-I collinear frequency-degenerate phasematching. As the pump, we used the 3rd 
harmonic of a Nd:YAG laser with the wavelength X p = 354.7 nm, pulse duration 18 ps, repe- 
tition rate 1 kHz and the energy per pulse up to 0.1 mJ. The coherence time of the pump was 
less than the pulse duration and amounted to 5 ps. To increase the parametric gain, the pump 
was focused into the crystal by means of a lens (LI) with the focal length 1 m. This resulted in 
the maximal peak intensity 10 GW/cm 2 . The gain value could be varied by changing the pump 
power, which was done by rotating a half- wave plate (HWP) in front of a Glan prism (GP). Af- 
ter the crystal, the pump radiation was reflected by a dichroic mirror (DM), which transmitted 
nearly all the down-converted light (97%), and its power was measured with a power-meter. 
The PDC radiation was filtered in the transverse wavevector (angle of emission) by means of 
an aperture (A) placed in the far-field zone. The emission angle selected this way was 0.17° 
(inside the crystal). All radiation passing through the aperture was focused, using a lens (L2), 
on the input slit of the HORIBA Jobin Yvon Micro HR monochromator equipped with a S YG- 
NATURE PDA CCD array. The spectral measurements were performed with the best possible 
resolution, which was 0.2 nm. 

In all our measurements, the spectral properties of PDC were studied depending on the para- 
metric gain values. In its turn, the parametric gain was found by fitting the dependence of the 
PDC output intensity I on the input pump power P (inset in Fig. []} by the function 

/ = /osinh 2 (K-v / P), (1) 

where To, k were the fitting parameters fT9l . The gain was measured for collinear frequency- 
degenerate regime. It is important that for such a measurement, the PDC radiation should be 
very well filtered both in the angle and in frequency, so that less than a single mode is registered. 
After the fitting, each pump power is put into correspondence to a certain value of the parametric 
gain G = K^/P. 

The dependence of the PDC spectral width (FWHM) on the parametric gain was measured 
for the range of gain values 3.9<G<6.5.At lower gain values, the PDC signal was too low to 
make spectral measurements while higher gain values could not be achieved with the available 
pump powers. The measured dependence is shown in Fig. O One can see that in the chosen 
range of gain values, the spectral width changes by 12%. Compared to the spectral width of 
low-gain PDC, the increase is 27%. Note that the effect of the gain on the spectral width is 
even stronger in the case of type-II or frequency-nondegenerate phase matching IfTTIl . The in- 
sets show the shapes of the spectra recorded at the lowest and highest gain values. The shapes 
are well fitted with the theoretical dependencies calculated using the approach of Ref. fT3l . 




Fig. 1. (color online) The experimental setup. PDC is generated in a type-I BBO crystal 
with the 3rd harmonic of Nd: YAG laser as a pump. The pump power can be changed using 
a halfwave plate HWP and a Glan prism GP. The pump is eliminated with a dichroic mirror 
DM. For PDC radiation, a narrow angular spectrum is selected with an aperture A. The 
spectral properties are studied using a spectrometer and a CCD camera. The inset shows 
the dependence of the PDC output intensity on the pump power. 



The same approach was used in calculating the theoretical dependence of the FWHM on the 
parametric gain (dashed line). One can observe perfect agreement between the theory and ex- 
periment. The only fitting parameter used in the theoretical dependence was the orientation of 
the crystal, which turned out to differ by 0.0025° from the one for exact collinear degenerate 
phase matching. This led to about 2 nm shift for the whole dependence, which was difficult to 
notice in the experiment. 

Here, it is worth mentioning that the 'raw' spectra obtained at the output of the spectrome- 
ter were always asymmetric, the long-wavelength part having smaller intensity. A typical dis- 
tribution is shown in Fig. [3] by red circles. It is interesting that similar distributions have been 
reported in many papers on degenerate (broadband) PDC - see, for instance, Ref . iTTSl . 

By analyzing the detection scheme (standard for such measurements), which includes the 
selection of a certain small solid angle and a certain wavelength interval, we come to the con- 
clusion that the spectra should be multiplied by a factor scaling as the fourth degree of the 
wavelength, A 4 . The reason is that the total number of photons is the product of the number 
of photons per mode (wavelength-independent in the vicinity of the degenerate wavelength) 
and the number of modes collected. At the same time, the number of transverse and longitu- 
dinal modes within a given solid angle interval A£l and a given wavelength interval AA scales 
as 1/A 4 . This can be understood from the following considerations. The number of modes is 
given by the ratio of the detection volume and the coherence volume l20l . 



The coherence volume is V co h = S co hCt co h, where S co h is the transverse coherence area, c is the 
speed of light, and t co h is the coherence time. The detection volume does not depend on the 
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Fig. 2. (color online) The dependence of the PDC spectral width on the parametric gain. 
The insets show the shapes of the spectra at gain values 3.9 (left) and 6.5 (right). The 
dashed line is exact calculation using the theoretical approach of Ref. JT3l . 




Fig. 3. (color online) The measured type-I PDC spectrum before (red circles) and after 
(black squares) the correction by A 4 factor. The spectrum corresponds to the gain value 
G = 6.48 ±0.05. 
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Fig. 4. (color online) Normalized variance of the difference between the PDC intensities 
at a given wavelength 702.4 nm and scanned wavelengths. Different curves correspond 
to different gain values: 6.4 (lower curve, black), 6.5 (middle curve, red), and 6.8 (upper 
curve, blue). 



wavelength, while the coherence area and time are wavelength-dependent, 

where Ak± is the interval of transverse wavevector variation and Acq is the frequency bandwidth 
corresponding to AA. As a result, the number of collected modes scales as A -4 , and this factor 
should be taken into account in the spectra. Figure [3] shows the shape of the spectrum before 
(red circles) and after (black squares) the correction. The correction also took into account the 
fact that in our CCD camera, the wavelength interval selected by the one pixel slightly depended 
on the wavelength. This lead to a 10% variation of the wavelength resolution along the CCD 
camera. 

In order to study the spectral correlations in high-gain PDC, we used the technique of noise 
reduction measurement [8, 9, 10], in which the observable under study is the variance of the 
difference of intensities for two beams. Correlation between two beams leads to a decrease in 
the noise of the difference intensity. We have studied the difference-intensity noise between 
beams at various pairs of wavelengths in the PDC spectrum. For different gain values, 10000 — 
20000 frames were made, each one collecting the intensity from 500 — 1000 pulses, and then 
processed. First, one wavelength was fixed (by fixing the pixel of the CCD camera) and the 
other one was scanned. The variance was calculated for the difference between the signals 
from the chosen pixel and all remaining ones. The results are shown in Fig.H In all spectra, the 
fixed wavelength ('signal') was chosen to be X s = 702.4 nm, while the other one ('idler'), A/, 
was scanned within the range from 700 nm to 720 nm. 

One can see that each spectral distribution of the variance contains two 'dips' and one peak. 
The left-hand 'dip' is caused by the classical correlations existing within a single mode of 
PDC radiation. The dip goes down to zero because the variance is identically zero for the case 
where the reading of a pixel is subtracted from itself. This point is an artifact and should not 
be considered. Another artifact is the oscillating noise seen on some (but not all) spectra in 



the vicinity of the left-hand 'dip'. The right-hand dip is caused by the nonclassical correlation 
between light beams scattered into parametrically conjugated modes, with the wavelengths 
related as A" 1 + A- -1 = A" 1 . 

The peak in Fig. [4] occurs at exactly the degenerate wavelength 709.3 nm. It appears be- 
cause at the degenerate wavelength, PDC radiation is single-mode squeezed vacuum, some- 
times referred to as 'light with even photon numbers', and this kind of radiation has peculiar 
statistical properties. Namely, single-mode squeezed vacuum manifests super-bunching even 
at strong parametric gain |2H . Its second-order normalized Glauber's correlation function is 
g( 2 ) = 3 + 1 /N, where N is the mean photon number per mode, and in the high-gain regime 
g( 2 ) = 3. Therefore, it has increased fluctuations, which causes the peak in Fig. 01 

The dips and the peak in Fig.|4]can be explained in more detail by recalling that the variance 
of the photon-number difference between the beams labeled s and i is l22ll 

Var(AT s - Nj) = (g%> - l)(N s ) 2 + {gf - 1 ) (N t ) 2 - 2{g<® -1)(N S ) (N f ) + (N s + N i ), (4) 

(2) (2) 

where g ss , g u are second-order normalized correlation functions for the signal and idler beams 

and gf^ is their cross-correlation function. The last term in d4j) represents the shot noise and 
is negligible in the high-gain regime. Note that Eq. ©, initially written for a single radiation 
mode, will not change if we assume that m independent modes are registered: the variance in 
the left-hand side, as well as the mean values in the right-hand side, will acquire a factor of 
m, and the factors containing the correlation functions in the right-hand side will change as 

g (2) _ i _> {g (2) _ 1)/m (Xgj. 

From the second term we see that the variance of the difference signal should have a peak if 
Aj corresponds to the degeneracy. Indeed, if (N s ) = (Ni) = N, then outside of the peak, where 
s\? = SsP — 2, the variance of the difference photon number is 2N 2 while in the peak, g^ = 3, 
and the variance is 3N 2 . In the 'autocorrelation dip' , A* ~ X s , g^ = g$ = 2, and Yar(N s —Ni) « 
2N, i.e., the variance of the difference photon number should be suppressed down to the shot- 
noise level. In the 'cross-correlation dip', where X s and A ; are conjug ated, gf) = 2 + l/N,and 
in theory, Var(A^ — Ni) should be zero. However, this complete suppression of noise can be 
observed only provided that three conditions are satisfied |[T0l : (1) the quantum efficiency is 
high; (2) a large number of transverse and longitudinal modes are collected and (3) the signal 
and idler modes are properly matched. In our experiment, the quantum efficiency is not more 
than 20%, only a few transverse modes are selected, and no proper mode matching is performed. 

We see that with the increase of the gain, the 'dips' in Fig.|4]become more pronounced. This 
is because the background is caused by the excess noise in the PDC intensity (the first two 
terms in ©), which scales quadratically with the mean photon number. At the same time, the 
minimum in the 'dip' is determined by the shot noise (last term in ©), scaling linearly with 
the mean photon number. As the gain grows, the mean photon number increases and the 'dip' 
becomes more visible. 

Another observation one can make from Fig.|4]is that the peak width is considerably less than 
the width of the 'dips'. Accurate fitting of all of them by Gaussian functions yields, in the case 
of the middle dependence in Fig. SI the values of 0.26 ± 0.06 nm, 0.45 ± 0.06 nm and 0.52 ± 
0.06 nm for the peak, the 'autocorrelation dip' and the 'cross-correlation dip', respectively. 
The widths of the dips coincide with the widths of the second-order correlation functions and 
indicate the spectra of frequency correlations. Theoretically, they should be given by the pump 
spectral width, which is Aco p = 1.25 THz. This corresponds to AA = 0.28 nm at wavelengths 
around 709 nm, but at high gain should be larger due to the nonlinear dependence of the PDC 
intensity on the pump intensity. 

As to the peak, its width, theoretically, should be narrower than the pump width by a factor 
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Fig. 5. (color online) Two-dimensional distribution of the covariance for the PDC intensities 
measured at signal and idler wavelengths (left) and its zoomed central part (right). The 
covariance is normalized to its maximum value. 



of 2. The observed widths of 0.26 ± 0.06 nm for the peak and 0.52 ± 0.06 nm for the dip are in 
agreement with this figure. 

More complete information about the frequency correlations in the spectrum of high-gain 
PDC can be obtained by measuring a two-dimensional distribution of some parameter charac- 
terizing the correlations, as a function of both wavelengths Ay, A;. As such a parameter, one can 
choose the difference- signal variance given by Eq. © but it is more convenient to measure the 
covariance of the signals at two various wavelengths, 

Cov(A^) = (N s Ni) - (N s )(Ni). (5) 

Note that the same strategy was used in Ref. fl5l , where the correlation coefficient (normalized 
covariance) was measured to characterize intensity correlations in the spectra of bright beams. 

The covariance is equal, up to a constant factor, to the third term in Q, since (gff — 
l)(N s )(Ni) = Cov(N s , Ni). It gives a full account of the pairwise correlations, which now mani- 
fest themselves as peaks (Fig. [5]). The 'autocorrelation' peaks occur on the straight line X s = A* 
and correspond to the classical correlations within a single PDC mode. The 'cross-correlation' 
peaks correspond to the signal-idler correlations and are on the line given by A" 1 + A- -1 = A" 1 . 
The distribution of these peaks is similar to the 'joint intensity' observed for two-photon 
light (161 [T71 HH : its cross-section (conditional distribution) shows the spectral range of corre- 
lations while its projections on the axes X s ,Xi (marginal distributions) show the PDC spectrum. 
The increased value of the correlation function for the collinear frequency-degenerate PDC also 
reveals itself in the vicinity of the point X s = X[ = 709.3 nm. 

In conclusion, we have performed an experiment aimed at the study of the spectral properties 
of high-gain PDC. In agreement with the theory, we observed the broadening of the spectrum 
with the increase in the parametric gain. We have shown that the dependence of the coherence 
volume on the wavelength causes the asymmetry of the measured PDC spectra, which can 
be removed by multiplying the spectra by a factor of A 4 . By measuring the variance of the 
intensity difference for two wavelengths, we saw the correlations in the PDC spectrum: both 
auto-correlations, due to the thermal statistics of each mode of the PDC radiation, and cross- 
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correlations, due to the pairwise statistics of signal and idler modes. The technique of noise 
reduction measurement allows us to resolve the spectral width of the intensity correlations. 
The measured variance distributions also show the super-bunching of PDC radiation at exactly 
degenerate wavelength. Finally, we have measured the distribution of the co variance for the 
PDC intensities at two wavelengths, both scanned within the range 650 — 760 nm. The obtained 
two-dimensional distributions show both the intensity correlations within a single mode of PDC 
radiation and the cross-correlation between the intensities at the signal and idler wavelengths. 
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